Negative Echo in the Density Evolution of Ultracold Fermionic Gases 
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We predict a nonequilibrium critical phenomenon in the space-time density evolution of a fermionic 
gas above the temperature of transition into the superfluid phase. On the BCS side of the BEC-BCS 
crossover, the evolution of a localized density disturbance exhibits a negative echo at the point of 
the initial inhomogeneity. Approaching the BEC side, this effect competes with the slow spreading 
of the density of bosonic molecules. However, even here the echo dominates for large enough times. 
This effect may be used as an experimental tool to locate the position of the transition. 
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Superconducting [or superfluid] fluctuations (SF) have 
a profound effect on the DC transport and the thermo- 
dynamic properties of superconductors, leading to the fa- 
mous Aslamazov-Larkin and Maki-Thompson (MT) cor- 
rections to the conductivity [1, 2]. Valuable informa- 
tion about the role of SF in the emerging correlations 
of the many-body wave functions is provided by the real 
time dynamics. Unfortunately, the latter can hardly be 
probed in solid state experiments. 

On the side of atomic physics, recent studies of 
fermionic ultracold atomic gases (UCAG) with inter- 
particle interaction that can be fine-tuned by exploit- 
ing Fano-Feshbach resonances (FFr)[3-6], have revealed 
pronounced effects and found clear signatures of super- 
fluidity [7-14] in the regime of the so-called BEC-BCS 
crossover [15]. In contrast to solid state systems, where 
it is hard to observe a dynamical response to external 
perturbations, in trapped UCAG it is possible to probe 
the actual evolution of a disturbance in real time and 
space. This offers, among other things, the exciting pos- 
sibility of discovering novel manifestations of the peculiar 
physics of SF. 

In this Letter, we analyze the nonequilibrium space 
and time evolution of density in the BEC-BCS crossover, 
and find that even above the critical temperature SF 
qualitatively affect the dynamics. The quantum kinetic 
equations (QKE) we derive here, allow us to analyze dy- 
namical effects in the critical region. 

Suppose that originally (at t = 0) the density is de- 
pleted in a small region inside the system. The time 
evolution of this density dip, illustrated by Fig. 1, turns 
out to be determined by how close the temperature T is 
to its critical value T c , and by the position of the sys- 
tem along the crossover. In the BCS limit, where T c 
is low, we predict the development of a sharp peak in 
the density of fermionic atoms at the origin (see Fig. 
la). The shape of the peak is critical: as T — > T c + the 
width shrinks as r 1 / 2 (r = (T — T c )/T c ). As the reso- 
nance detuning is further lowered, T c increases and cor- 
rections to the pure BCS behavior become important: 
the bosonic molecules depleted initially produce a dip 
that dominates the fermionic peak for small times and 
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FIG. 1: (Color online) Density profile 6pp(x.,t), Eq. (16), in 
arbitrar y units as a function of the scaled distance from the 
origin *J irb/y(E-i_)mr/t, for values of r shown on the curves. 
The inset sketches the total density evolution in regions III- 
IV (a) and II (b). However, due to the exponential decay of 
the bosonic dip, even in region II one always ends up with an 
algebraically decaying fermionic peak (a-3) for large times. 



r's (see Fig. lb and region II of Fig. 2). Eventually, how- 
ever, the bosons decay and the fermionic peak emerges. 
We will find that tuning the system towards the bosonic 
side of the crossover reduces the interval of r's where the 
fermionic density peak is observable. Eventually, in the 
BEC regime only the bosonic dip remains. To measure 
these effects, we propose the experiment shown in Fig. la 
with fermionic UCAG close to the superfluid transition, 
following the experimental procedure introduced by [16], 
e.g. by applying the repulsive optical dipole force of a 
strongly focused off-resonant laser beam. We stress that 
the negative echo effect is not limited to the far BCS 
limit and is observable in the region that is accessible to 
current experiments. 

The profile of the density peak formed by fermions in 
the BCS side of the crossover, is plotted in Fig. 1. The 
microscopic mechanism at the root of this phenomenon is 
the Andreev reflection of fermionic atoms off the fluctua- 
tions of the superfluid order parameter. The wavelength 
of such fluctuations diverges at T — * T c + . It is well known 
[17] that an electronic excitation, with sub-gap energy, in- 
cident on a superconductor/normal metal interface from 
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the normal side, will be reflected back as a hole, almost 
opposite to its original direction of motion, creating a 
Cooper pair in the superconductor. Although the su- 
perconducting gap vanishes above the transition point, 
close to T c there appear critically large fluctuating su- 
pcrfluid islands which transform incident nonequilibrium 
particle-like excitations into holes, and vice versa. Thus 
the time evolution shown in Fig. la can be understood 
as follows. The initial density depletion can be thought 
of as an excess of holes, which at t = start to move 
radially outwards from the origin, while the fluctuating 
condensate reflects them back towards the origin, in the 
form of atoms rather than holes. Possible deviation of 
this reflection from exact backscattering is determined 
by the typical size of a condensate island, which diverges 
as t -1 / 2 . As a result fcrmions accumulate at the ori- 
gin forming a critically sharp density peak, the width of 
which scales as t 1 / 2 (see Eq. (16)). 

On the other hand, the density dip of Fig. lb, is due 
to the slow propagation of the bosonic component of the 
initial nonequilibrium disturbance away from the origin 
(see Eq. (17)). 

The fermionic atoms interacting via a FFr are de- 
scribed by [18-20] 
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e a =p 2 /(2m,) + U a (r,t); e b = p 2 /(4m)+U b {r,t), (lb) 

where a^(r), [d a (r)] create [annihilate] fermionic atoms 
with mass m and spin a = ±1/2, while ft(r), [b(r)} create 
[annihilate] the bosonic molecules. U a ^) 1S the sum of the 
confining and perturbing potentials acting on an atom 
(molecule). 

Given the energy of relative motion of two atoms, e, 
their s-channel scattering amplitude can be written as 



/o(e) = h(me) 1/2 7(e) [e - e - h(e)]~ 
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with the resonance width 7(e) 



and physical 



position e . Physical results expressed in terms of e and 
7 are well defined [21]. 

QKE for the normal phase are derived (see e.g. 
Ref. [22]) by finding a semiclassical approximation for 
the time ordered Green functions. We combine the lat- 
ter into 2x2 matrices, 
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and G K I R I A are 2x2 matrices in the spin space. 

We limit ourselves to the condition of narrow reso- 
nance, (see Ref. [23] for the experimental realization) 

7(E) <C £, E = max (T, e F ) , 2me F = [3ir 2 h 3 p a ] f , (4) 



where p a is the atomic density in the absence of the reso- 
nance, and T is the temperature [24]. For r S> Gi, where 
Gi <C 1 is the Ginzburg parameter [25], condition (4) 
justifies the one loop approximation 
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for the derivation of the QKE. Here the structure of the 
vertices in the Keldysh space (3) is defined by dkim = 
Sk2$im + 8126km + 5 m2 8ki - 2Si 2 S k2 S m 2. As usual [22], the 
derivation of the QKE requires obtaining the equation 
for the Wigner transforms for the Green function (1-1.2 = 
r±5r/2, t ia = t±5t/2): 
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and keeping leading terms in h expansion. [Hereinafter, 
we use the notation X = (p, x), x = (r, t).] 

We parametrize the Keldysh components of matrices 
(3) in terms of the distribution functions n(X),N(X): 
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with h(X) a 2 x 2 matrix in the spin space. Hereinafter, 
[-,•]_!_ means anticommutator. Using Eqs. (5)-(7), we de- 
rive the QKE (details will be published elsewhere): 



d t h+- 
h 



d t N + {n,N} p = - 



dpi 

P J {2nh} 
dp 1 



Si(p,Pi;x); (8a) 



2 (2irh)- 



TvSt (p- Pl , Pl ;x),(8b) 



where [•,•] stand for the commutator and the Pois- 
son brackets for arbitrary matrices A, B are defined as 
2{A,B} p = [V p A,V r B] + - [V r A,V p B] + . Equation 
(8b) is valid when the bosonic spectral width is smaller 
than temperature, and thus is not valid in the far BCS re- 
gion (region IV of Fig. 2). However, it can be shown that 
in this region it suffices to use the equilibrium bosonic 
distribution (sec Eq. (13)) in Eq. (8a), and Eq. (8b) is 
not needed. 

In Eq. (8b), the energies of the fermionic and bosonic 
excitations, £ and il, depend on the momentum and the 
coordinate through the bare energies (lb) and also on 
the distribution functions of the excitations themselves. 
Due to a possible non-equilibrium spin population, the 
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spectrum £ may posess a spin space structure: 
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where notation A = fl(p 1 + p 2 ,x) — £{P\, X ) — £{p 2 , x ) 
is introduced and overline signifies the time reversal op- 
eration in the spin space: A = a y A T a y . The spectrum 
of the bosonic excitations acquires a shift 
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FIG. 2: Relevant regions of parameter space. Region I is the 
superfluid phase. The QKE (8b) is valid in II and III. Our 
results are valid above and close to the critical line (black 
curve). In III and IV the MT correction, Eq. (16), dominates 
at all times, whereas in II the bosonic dip, Eq. (17), dominates 
for t < t, . 



which for a sharp fermionic distribution function, di- 
verges logarithmically as in conventional BCS. 

Real collision processes defining the irreversible evolu- 
tion in Eq. (8b) are described by (we omitted x in the 
argument of the distribution functions) 
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T = N(p 1 +p 2 ) 1 - n(p x ) -n(p ; 
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The QKE scheme (7)-(10) is justifiable as long 
as the temporal [spatial] variation of the distribu- 
tion functions is smooth on quantum scale, h/T 
[7i/vWTmax(l, yJep/T)} [24]. Under these conditions, 
the expressions for the collision integral and spectral 
renormalization are local in time and space [26] . 

For densities of fermions, p a {x), bosons pb{x), and 
spin, S(x), we have 

[p a ,2S,p b ] (x) = Jj0-[Trh,Tr*n,N] (X). (11) 

Quantum Kinetic Equations (8) respect the symme- 
tries of model (1): conservation of mass, energy, spin 
rotation invariance, and when VC/ a- 6 = 0, Galilean in- 
variance. The latter is implemented by 



n(p, r) — > h(p — mv, r — vt); 
N(p, r) — > N(p — 2mv, r — vt); 
£(p, r) — > £(p — mv, r — vt) + v ■ p — mv 2 /2; 
f2(p, r) — > VL{p — 2mv, r — vt) + v ■ p — mv 2 . 
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Furthermore, for dtU a .b = the equilibrium distributions 
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where T{x) = l/(e x + 1); B{x) = l/{e x - 1), solve 
Eqs. (8b), with arbitrary constants, fJ,,T, and an arbi- 
trary traceless Hermitian 2x2 matrix, h. The latter 
describes a possible stationary imbalanced spin popula- 
tion, which will be discussed elsewhere. 

Transition temperature, T c , is given by [27-31] 
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In the BCS regime, e — 2ep ^> 7(ej?), this yields T c cx 



e F e 2 "'< 2£ f> , and in the BEC regime, (e < 0), T c w 
Tbec — 0.218ejr. The equilibrium bosonic spectrum, 
in the interval T c Gi < O(q) - 2/i < T [25], is ft(q) = 
2[i + c + bq 2 , with 
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Here, £(3) ~ 1.202 (Riemann (- function) , and 
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In the BCS limit T c < ^, E 1<2 -> 2^ w 2e F . 

We now consider the evolution of a small density dis- 
turbance as shown in Fig. 1. A realistic experimental 
procedure [16] for creating such a disturbance is by ap- 
plying a potential C/ Q (x) > 0, much narrower than the 
size of the trap, for a time long enough for the system 
to equilibrate, with n(X,t = 0) = T f §lskiI^EsM \ 

At time t — 0, the potential is removed abruptly, leav- 
ing a small density disturbance. Next, we linearize 
the QKE and the inital distributions in the deviation 
from equilibrium. We solve the linearized fermionic ki- 
netic equation (8a) perturbatively in the collision inte- 
gral, Eq. (10), but treat the equilibrium decay and spec- 
tral renormalization of the bosons non-perturbatively, 
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consistent with our one-loop approximation. To zeroth 
order, the nonequilibrium part of the fermionic distri- 
bution evolves as that of free particles: 5n^(X,t) = 
5n(p, x—pt/m, t = 0), giving rise to an expanding spher- 
ical density wave. Similarly for bosons, except for their 
exponential decay: 5N^{X,t) = SN(q,x - 2bqt 1 t = 



0) exp 



( C +bqf 2 ) 7 (2 M ) ■ 
2T C u 



Deep in the BCS regime (regions 

III and IV in Fig. 2) the nonequilibrium bosons decay fast 
and can be neglected. Iterating the fermionic QKE, we 
find that in this regime, the correction to the density close 
to the origin is dominated by the same terms in the col- 
lision integral (10) that give rise to the MT correction to 
the conductivity of conventional superconductors. The 
latter can also be understood as the effect of Andreev 
reflection of electrons by superconducting fluctuations. 
For times t 3> tp = T 2*V ; the density peak normalized 



by the total number of initially removed fermions is de- 
scribed (for r«W^) by 
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In the BCS limit, this result is valid even for a broad 
resonance, i.e. not restricted by the condition Eq. (4). 

In contrast to the exponential decay of any nonequi- 
librium deviation in the bosonic distribution, this sharp 
peak in the fermionic density decays only algebraically 
(one can show that this conclusion is not an artifact of 
perturbation theory), so that for large times the lat- 
ter always dominates. For times smaller than t* = 
7 (2 M )^(g 1 ) ' an d close enough to T c , the increasingly 
slow moving nonequilbrium bosons removed by the ex- 
ternal potential [32] produce a density dip at the origin, 
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where v 2 = 8lry(Ei)/ir. This formula is valid for t 
ts = h/ \J 'j(Ei)T c t and r <C t\fbT c , and dominates the 
fermionic peak Eq. (16), reversing the sign of the density 
deviation at origin (see Fig. lb). However, in the BCS 
regime, t* is suppressed, and the peak described by for- 
mula (16) dominates at all times. The region where MT 
correction, Eq. (16), dominates is sketched in Fig. 2. 

The critical behavior that was found in the response 
of the system to an external force may be of valuable 
use in current experiments on UCAG. It can provide for 
instance a new experimental tool to determine the tran- 
sition temperature in the BCS side (where the formation 
of vortices [12] is not a strong enough mechanism), by 
applying a small density disturbance on the system and 
measuring nondestructively its subsequent evolution. 

In summary, we undertook the study of the real time 
density evolution of a composite fermion-boson system in 
the BEC-BCS crossover, close to the transition point in 



the normal phase. We found that in the far BCS regime 
a localized depletion of density will evolve into a sharp 
diverging peak where it was created, as a manifestation 
of critical fluctuations of the supcrfluid condensate. This 
phenomenon can be exploited in the experimental study 
of the crossover in fermionic UCAG's. 

We are grateful to G.V. Shlyapnikov for inspiring con- 
versations that motivated this work, and to G. Catclani, 
A. Andreev, and L. Glazman for useful remarks. 
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